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We have established detailed magnetoelectric phase diagrams of 
(Eu .595Y .40 5 )i-xTb ;c MnO 3 (0 < x < 1) and (Eu,Y) 1 _ x Gd ;c Mn03 (0 < x < 0.69), 
whose average ionic radii of i?-site (R: rare earth) cations are equal to that of Tb 3+ , in order 
to reveal the effect of rare earth 4/ magnetic moments on the magnetoelectric properties. 
In spite of the same i?-site ionic radii, the magnetoelectric properties of the two systems 
are remarkably different from each other. A small amount of Tb substitution on R sites 
(x ~ 0.2) totally destroys ferroelectric polarization along the a axis (P a ), and an increase 
in Tb concentration stabilizes the P c phase. On the other hand, Gd substitution (x ~ 0.2) 
extinguishes the P c phase, and slightly suppresses the P a phase. These results demonstrate 
that the magnetoelectric properties of i?Mn03 strongly depend on the characteristics of the 
rare earth 4/ moments. 
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In recent years, a new class of multiferroics, which have both magnetic and ferroelec- 
tric orders, has been attracting increasing interest. i?MnC>3 (R: rare earth) of orthorhombic 
perovskite structure is one of the most typical examples. In R = La - Eu, i?Mn03 has the 
^4-type antiferromagnetic ground state of Mn 3d spins. 1 ) As the perovskite i?-site ionic ra- 
dius decreases, the magnetic ground state changes from A-type antiferromagnetic order to 
noncollinear transverse spiral antiferromagnetic order around R = Gd, owing to magnetic 
frustration between the nearest neighbor (NN) ferromagnetic and next NN antiferromagnetic 
interactions among Mn spins. 2 ~ 6 ) In i?Mn03 where R = Tb or Dy, noncollinear transverse spi- 
ral antiferromagnetic order breaks the inversion symmetry through the inverse Dzyaloshinskii- 
Moriya (DM) interaction, resulting in the appearance of ferroelectricity. 3 ' 5 ' 7 ' 8 * 1 The spin he- 
licity or vector spin chirality of the spiral spin structure can be changed by external magnetic 
fields, then the direction of the ferroelectric polarization due to the spin chirality can also 
be controlled by the fields. 9 ) As the .R-site ionic radius decreases from R = Dy to R = Ho - 
Lu, and Y, the ground state changes from noncollinear spiral antiferromagnetic to collinear 
.E-type antiferromagnetic order. 1 ' 10 ) i?Mn03 with E-type antiferromagnetic order also ex- 
hibits ferroelectric behavior that is caused not by the inverse DM interaction but by inverse 
Goodenough-Kanamori interaction, i.e., the symmetric exchange striction. 10 ' u ) 

The magnetoelectricity of i?Mn03 is largely affected not just by the magnetic frustration 
of Mn 3d spins but also by 4/ magnetic moments of i? 3+ . 12 ~ 17 ) In our previous works, 18 ) to 
clarify the effect of 4/ moments, we made a comparative study of TbMn03, Yo.3iGdo.6gMn03, 
and Euo.sgsYo^osMnOs, whose average ionic radii of R sites are equal to that of Tb 3+ . Al- 
though they have the same ionic radii of -R-site cations, the magnetoelectric properties of 
these three compounds are quite different. In 4/-moment-free Euo.595Yo.405Mn03, 13,19 ) fer- 
roelectric polarization along the c-axis (P c ) emerges below TpEc = 25 K; P c disappears at 
TpEa = 23 K, below which ferroelectric polarization is observed parallel to the a axis (P a ) 
[Fig. 1(a)]. In TbMnOs, the polarization flop does not occur without a magnetic field, and the 
polarization of the ferroelectric ground state is along the c-axis. 2,8,9 ) The difference between 
their magnetoelectric properties is attributed to the presence of Tb 3+ 4/ moments. Tb 3+ 
has an anisotropic 4f moment, because of its large orbital angular momentum (L = 3). The 
ordered Tb 3+ moments below 7 K are along either of two Ising axes within the ab plane at 
57° off the 6-axis. 20 ) Thus, in TbMnOs, magnetic coupling between 4/ Ising-like moments and 
Mn 3d spins has an effect similar to that of applying external magnetic fields to the Mn 3+ 
sublattice. Hereafter, we refer to these effective magnetic fields as "internal magnetic fields" . 
The 4/-moment effect of Gd 3+ is opposite to that of Tb 3+ . In Y .3iGd .69MnO 3 , only the P a 
phase is found in a zero magnetic field. 14 ' 18 ) The 4/ moment of Gd 3+ is of the Heisenberg type 
without anisotropy, because L = 0. Consequently, the Mn 3+ sublattice of Yo.3iGdo.6gMn03 
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does not feel the internal magnetic field along the a-axis, which develops P c . 

As mentioned above, the rare earth 4/ moment has a significant effect on the magne- 
toelectric properties of i?MnC>3, but the effect is not well understood yet. In this study, 
to reveal the contribution of 4/ moments to the magnetoelectric properties of i?Mn03, we 
have established detailed magnetoelectric phase diagrams of (EU0.595Y0. 405)1-3; Tb^MnG^ and 
(Eu,Y)i_ x Gd x Mn03, whose average i?-site ionic radii are fixed at that of Tb 3+ . 

2. Experiment 

We prepared a series of (Euo.595Yo.405)i-xTb x Mn03 (0 < x < 1) and (Eu,Y)i_a;Gd x Mn03 
(0 < x < 0.69) crystals with several substitution levels, where the average ionic radii of R- 
site are equal to that of Tb 3+ . In the latter system, therefore, the ratio of Eu to Y is tuned 
appropriately depending on Gd concentration x. It should be noted that x = 0.69 is the limit 
of Gd concentration for retaining the average i?-site ionic radii at that of Tb 3+ . The average 
ionic radii were calculated on the basis of the Shannon's ionic radius table. 21 ) The starting 
materials, EU2O3, Y2O3, Tb4C>7, Gd203, and Mn304 were weighted to a prescribed ratio and 
mixed by using ethanol. The mixture was heated at 1073 K in air with a few intermediate 
grindings, and sintered at 1473 K in air. The resultant powder was formed to a cylindrical 
shape (6 mm (f> x 100 mm) with use of a hydrostatic pressure of 300 ~ 350 kgf/cm 2 to 
make a feed rod. Single crystalline samples were grown under Ar atmosphere (~ 3 atm) at a 
feed rate of 3 ~ 6 mm/h by a floating zone method. We performed X-ray diffraction (XRD) 
measurements on the obtained crystals at room temperature. The Rietveld analysis of powder 
XRD data revealed that all the samples have an orthorhombic Pbnm structure without any 
impurity phases or any phase segregation. All the single crystalline samples used in this study 
were cut along the crystallographic principal axes into a rectangular shape using an X-ray 
back-reflection Laue technique. The dielectric constant and the ferroelectric polarization in 
magnetic fields were measured using a temperature-controllable cryostat equipped with a 
superconducting magnet that provides a magnetic field up to 8 T. The dielectric constant 
measurement was performed with an LCR meter at a frequency of 10 kHz (Agilent, 4284A). 
The pyroelectric current to obtain the spontaneous electric polarization was measured in a 
warming process at a rate of 4 K/min after the samples were cooled from 60 to 5 K in a poling 
electric field of 300 ~ 500 kV/m. The magnetic properties were measured using a commercial 
apparatus (Quantum Design, PPMS-9T). 

3. Results 

Let us start with the 4/-moment-free compound, Euo.595Yo.405Mn03. 13 ' 16 ) Figure 1 shows 
the temperature dependence of (a) the ferroelectric polarization in a zero magnetic field and 
(b) the magnetization in H = 0.5 T. P c emerges below Tfe c = 25 K, and the magnetization 
parallel to the c-axis (M c ) shows a slight decrease at the same time. However, no anomaly is 
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clearly discerned in the magnetization parallel to the a-axis (M a ) or the o-axis (Mb) at Tfec- 
The P c phase is replaced by P a . That is, the direction of the ferroelectric polarization changes 
from the c-axis to the a-axis at Tfeo- This polarization flop is caused by the magnetic phase 
transition from be spiral (P c ) to ab spiral (P a ) antiferromagnetic order, which accompanies a 
decrease in M a and an increase in M c . The magnitude of P a is much larger than that of P c . 
This is because the polarization flop from the c-axis to the a-axis immediately below Tfec 
before the saturation of P c . The decrease in Mb below 47 K is due to an incommensurate 
antiferromagnetic ordering with a collinear sinusoidal modulation along the 6-axis, which is 
considered unrelated to the magnetoelectric properties. 

Tb doping into Euo.595Yo.405Mn03 greatly affects the magnetoelectric properties. In 
(Euo. 595 Yo. 405)0. 9Tbo.iMn03, TpEa shifts to lower temperatures by 8 K, and consequently 
the P c region is broadened compared to that of Euo.595Yo.4osMn03 [Fig. 1(c)]. M c shows a 
clear decrease at Tfec (= 25 K), and then increases below TpEa (= 15 K) [inset of Fig. 1(d)]. 
These changes are also consistent with the magnetic transition from the be spiral (P c ) to 
the ab spiral (P a ). Note that the magnetic anisotropy of (Euo. 595 Yo. 405)1-1 Tb x Mn03 is much 
larger than that of Euo.sgsYo^osMnOs. The observed M a is larger than Mb and M c over the 
whole temperature region, and magnetic anisotropy develops with decreasing temperature or 
increasing Tb concentration. Figure 2 shows the temperature dependence of the ferroelectric 
polarization of (Euo. 595 Yo. 405)1-0; Tba;Mn03 with < x < 0.2 in a zero magnetic field. With 
increasing Tb concentration, the P c phase develops abruptly, and Tfec rises slightly. In con- 
trast, the P a phase is more fragile against Tb doping, and it totally disappears at a Tb-doping 
level of x = 0.2. In x > 0.2, (Euo.595Yo.405)i-a;Tb x Mn03 shows magnetoelectric behavior sim- 
ilar to that of TbMnC>3. As we previously reported, this is probably because of the internal 
magnetic field arising from magnetic coupling between Mn 3d spins and Tb 4/ moments. We 
summarize these results in the magnetoelectric phase diagram shown in Fig. 3(a). 

The effect of Gd substitution contrasts sharply with that of Tb substitution. In 
(Eu,Y)o.gGdo.iMn03, the P c phase is squeezed into a narrow region of T = 23 K - 24 K, 
and the magnetic anisotropy is quite small compared to that of (Euo.595Yo.4os)o.9Tbo.iMn03 
[Figs. 1(e) and 1(f)]. Figure 4 shows the temperature dependence of P a and P c in 
(Eu,Y)i_ :E Gd :E Mn03 with < x < 0.2 in a zero magnetic field. As seen from Fig. 4, P a 
is relatively robust against Gd substitution, and it still remains up to x = 0.69. However, 
TpEa shifts to lower temperatures with increasing Gd concentration. In the end material, 
Yo.3iGdo.69Mn03, TpEa is reduced by 8 K. On the other hand, P c is not very robust against 
Gd-doping, and it is totally extinguished at x > 0.2. The magnetoelectric phase diagram of 
(Eu,Y)i_ x Gd x Mn0 3 with < x < 0.69 is displayed in Fig. 5. 
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First, let us consider the effect of Tb substitution. As shown in the magnetoelectric phase 
diagram of (Euo.595Yo.4os)i-xTb :r Mn03 (0 < x < 1) [Fig. 3(a)], with increasing Tb concentra- 
tion, TpEa falls sharply, and the P a phase completely disappears at x = 0.2. On the other hand, 
the P c phase exists for all values of x, and Tfe c rises slightly with increasing Tb concentration. 
We display in Fig. 3(b) the magnetoelectric phase diagram of Euo.595Yo.405Mn03 as a function 
of applied magnetic field parallel to the a-axis (H a ) for comparison with the magnetoelectric 
phase diagram of (Eu . 5 95Y .405)i-xTb :r MnO3. 13) In Eu . 5 95Yo.405Mn0 3 , with increasing H a , 
TpEa decreases, and P a is extinguished in H a > 5 T. The P c phase exists over the whole 
H a region, and TpEc is almost invariable against H a . This is because the application of H a 
destabilizes the ab spiral that induces the P a phase. Namely, H a perpendicular to the be spiral 
plane induces a conical state and gains Zeeman energy. As clearly seen in Figs. 3(a) and 3(b), 
the magnetoelectric phase diagram of (Euo.595Yo.405)i- :r Tb r Mn03 is qualitatively similar to 
that of Euo.595Yo.405Mn03 as a function of H a . From this similarity, one can conclude that 
the effect that Tb substitution has on the Mn 3+ sublattice as similar to the effect the appli- 
cation of H a has. The enhancement of M a induced by Tb substitution strongly supports the 
above conclusion [Fig. 1(d)]. The P a phase of Euo.595Yo.405Mn0 3 is destroyed by an external 
magnetic field of H a = 5 T, while that of (Euo.595Yo.405)i- a; Tb a ,Mn03 disappears at x = 0.2. 
Thus, the internal magnetic field of (Euo.595Yo.405)o.sTbo.2Mn03 can be estimated to be 5 T. 

Next, we discuss the effect of Gd substitution. As Figs. 3 and 5 show, the magneto- 
electric phase diagram of (Eu,Y)i_ x Gd x Mn03 (0 < x < 0.69) is quite different from that of 
(Euo.595Yo.405)i-xTba;Mn03. With increasing Gd concentration, both TpEc and TpEa decrease, 
and the P c phase cannot be observed at x > 0.2. This difference originates from the difference 
in the characteristics of the 4/ moments of Gd 3+ and Tb 3+ , which greatly affect the magnetic 
anisotropy of i?Mn03. The magnetic anisotropy of (Eu,Y)i_ x Gd x Mn03 is much smaller than 
that of (Euo.595Yo.405)i-xTb x Mn03 [Figs. 1(d) and 1(f)], because of the isotropic 4f moment 
of Gd 3+ . The nearly isotropic magnetization of (Eu,Y)i_ x Gd x Mn03 indicates that magnetic 
coupling between Gd 4/ moments and Mn 3d spins does not enhance the magnetization in a 
particular direction, which is in contrast to the case of (Euo.595Yo.405)i- a; Tb a: Mn03. Therefore, 
it does not act as an internal magnetic field. 

Another remarkable difference between the two systems is that both ferroelectric 
transition temperatures (TpEa and Tfe c ) decrease in (Eu,Y)i_ rr Gd z Mn03 but not in 
(Euo.595Yo.405)i-zTb z Mn03. The decrease in the ferroelectric transition temperatures is prob- 
ably due to i?-site randomness arising from solid solution of the i?-site ions and/or magnetic 
randomness arising from random distribution of Gd 3+ or Tb 3+ in R sites. Generally, the R- 
site randomness can be measured by the variance a 2 = T,(gir 2 — rf^j, where gi, ri (i = 1 ~ 3), 
and vr are the fractional occupancy, the ionic radius of the i-th ion, and the average ionic 
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radius of R sites, respectively. 22 ' 23 ** is equal to the ionic radius of Tb 3+ . a 2 , i.e., the -R-site 
randomness of both systems, becomes smaller with increasing Gd or Tb concentration. If the 
i?-site randomness effect were dominant, Gd or Tb substitution could raise the ferroelectric 
transition temperatures. However, Gd substitution actually reduces these transition tempera- 
tures. Spiral antiferromagnetic order arises from subtle balance between the NN ferromagnetic 
and the next NN antiferromagnetic interactions so that it would be rather sensitive to mag- 
netic randomness. As a result, Gd substitution reduces TpEc an d ?feo- On the contrary, in 
the Tb-substituted sample, the internal magnetic field effect, which stabilizes the P c phase, 
is larger than magnetic randomness effect. Consequently, TpEc of the Tb-substituted sample 
rises slightly. 

5. Conclusion 

In this study, we have established detailed magnetoelectric phase diagrams of 
(Euo.595Yo.405)i-zTb z Mn03 and (Eu,Y)i_ x Gd x Mn03 to clarify the effect of rare earth 4f 
moments on the magnetoelectric properties; the average i?-site ionic radii of the two sys- 
tems equal that of Tb 3+ . The parent compound Euo.595Yo.405Mn03 without 4/ moments 
undergoes the P c transition at 25 K, and then the P a transition at 23 K. Tb substitution 
has the similar effect on the Mn 3+ sublattice as applying H a (we refer to this effective mag- 
netic field as the internal magnetic field). As a result, the magnetoelectric phase diagram of 
(Euo. 595 Yo. 405)1-3; Tb^MnOs is qualitatively similar to that of Euo.595Yo.405Mn03 as a function 
of H a . On the other hand, in (Eu,Y)i_a;Gd x Mn03, the internal magnetic field is ineffective 
because of an isotropic 4/ moment of Gd 3+ . With increasing Gd concentration, TpEa and 
TpEc are decreasing, and the P c phase is totally extinguished in x > 0.2. Probably, magnetic 
randomness coming from random distribution of Gd 3+ in i?-site is the main cause of the de- 
crease in the ferroelectric transition temperetures. The present results reveal the significant 
contribution of 4f moments to the magnetoelectric properties of -RMJ1O3, and in addition, 
provide significant information to control the magnetoelectric properties using rare earth 4/ 
moments. 
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Fig. 1. (Color online) Temperature dependence of the ferroelectric polarization along the a (P Q ) 
and c (P c ) axes (upper) and magnetization along the a (M a ), b (Mb), and c (M c ) axes 
(lower) of Euo.595Y .405Mn0 3 [(a) and (b)], 13 > (Eu .595Yo.405)o.9Tbo.iMn0 3 [(c) and (d)], and 
(Eu,Y)o.gGdo.iMn03 [(e) and (f)]. Arrows in (b) show ferromagnetic transition temperatures of 
EU0.595Y0.405M11O3. Inset of (d) shows magnification of M c in (Euo. 595 Yo. 405)0. 9Tbo.iMn03. Ap- 
plied magnetic field of H = 0.5 T was directed along each principal axes. 
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Fig. 2. (Color online) Temperature dependence of the ferroelectric polarization along the a (P a , up- 
per) and c (P c , lower) axes of (Euo.595Yo.405)i-a;Tb a ;Mn03 with < x < 0.2 in a zero mag- 
netic field. 
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Fig. 3. (Color online) Magnetoelectric phase diagram of (E110.595Y0. 405)1-2: Tb^MnOa with < x < 1 
(a), and that of Euo.595Y .405Mn03 as a function of applied magnetic field parallel to the a-axis 
(H a ) (b). 13 ^ Abbreviations denote the paramagnetic (PM), paraelectric (PE), antiferromagnetic 
(AFM), and ferroelectric (FE) phases. Data points are obtained from the measurements of the 
dielectric constant. 
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Fig. 4. (Color online) Temperature dependence of ferroelectric polarization along the a (P a , upper) 
and c (P c , lew 
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Fig. 5. (Color online) Magnetoelectric phase diagram of (Eu,Y)i_ :E Gd a ;Mn03 with < x < 0.69. 
Abbreviations denote the paramagnetic (PM), paraelectric (PE), antiferromagnetic (AFM), and 
ferroelectric (FE) phases. Data points are obtained from the measurements of the dielectric con- 
stant. 
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